1 



© 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



© Publication number: 



0 295 849 

A2 



© 


EUROPEAN PATENT 


APPLICATION 


© 


Application number: 88305383.7 


© Int. CI. 4 : G02C 7/02 


© 


Date of filing: 13.06.88 








Priority: 19.06.87 US 63971 


© 


Applicant: SOLA INTERNATIONAL HOLDINGS, 








LTD. 


© 


D^tfl nf ouhlieation of aDOlieation* 




P.O. Box 244 




21 12 88 Bulletin 88/51 




Lonsdale Morphett Vale South Australia 








5162(AU) 


Designated Contracting States: 


© 






DE FR GB IT 


Inventor: Barkan, Eric F. 








3426 Balboa Street 








San Francisco California 94121(US) 








Inventor: Sklar, David H. 








363 Sixth Avenue 








San Francisco California 94118(US) 






© 


Representative: Jenkins, Peter David et al 








Page White & Farrer 5 Plough Place New 








Fetter Lane 








London EC4A 1HY(GB) 



© Improved progressive lens and manufacture thereof. 

© A method for optimizing the performance of an ophthalmic progressive lens having a distance viewing zone, 
a near viewing zone, and an intermediate zone, and the resulting optimized lens are provided. The surface of the 
lens to be optimized, which has at least one quantitative performance characteristic, is represented by a 
mathematical base function which is combined with a selected optimization function. The optimization function 
includes at least one adjustable parameter and may have at least one boundary condition imposed on it along 
the boundary of a selected subregion of the lens surface. A quantitative merit figure of the combination of the 
base and optimization functions, which is related to the at least one performance characteristic by a selected 
merit function, is optimized through the at least one adjustable parameter of the optimization function. Through 
the optimization, there is obtained an improved surface function which represents the surface of the improved 
progressive lens. A lens surface responsive to the improved surface function is then formed. 
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IMPROVED PROGRESSIVE LENS AND MANUFACTURE THEREOF 

The present invention relates to a method for manufacturing an optimized opthalmic progressive lens, 
and in particular which has an improved optical performance compared to an existing progressive lens or 
progressive lens design. The performance is expressed as a quantitative merit figure and is improved by 
selectively varying the shape of the lens surface to optimize the resulting lens 1 merit figure. The present 

s invention also relates to a multi-focal progressive lens. 

The method of the present invention may be applied to an ophthalmic progressive lens having distance, 
intermediate and near viewing regions. Methods of making progressive lens surfaces are described in the 
Applicants' EP-A-01 32955 published 13th February 1985 and the disclosures of that publication are hereby 
incorporated by refereace in the present application as if fully set forth. Lenses made according to methods 

70 described in application Serial No. 6/516,366 have been used publicly and sold in the United States more 
than one year prior to the date of the present application. 

"Improving" or "optimizing" a lens design involves some criterion of merit Depending on the purpose 
intended for the lens, the merit criterion may weight various quantifiable lens characteristics more heavily 
than others. However, once such relative weightings have been chosen, a weighted combination of 

75 quantitative measures of lens performance can be formulated. An evaluation rule can thus be constructed 
that yields a single number which reflects the overall "quality" of the design according to the chosen 
weightings. Such a rule is often called a "merit function" and its value is often called a "merit figure". The 
merit function can be constructed such that either smaller or larger values of the merit figure indicate better 
lens performance, in this context, an "optimum" lens is one that minimizes or, if appropriate, maximizes the 

20 value of the merit function. A variety of such methods for minimizing or maximizing the values of such merit 
functions in general lens design are known in the prior art. The CODE V program available from Optical 
Research Associates, 550 North Rosemead Boulevard, Pasadena, California, is one such general lens 
design tool. The mathematics of such optimization procedures are discussed in, for example, Hamming, 
R.W., Numerical Methods for Scientists and Engineers, 2d ed. t Dover Publications, New York 1973, Chapter 

25 43. 

Important quantitative measures of lens performance are astigmatism, orthoscopy, and mean curvature. 
These are optical measures which may be computed in two somewhat different ways. The first way, which 
may be termed surface measure, uses the properties of an optical surface unmodified by any effects of ray 
obliquity that may occur in actual use of the surface. It is known to compute such surface measures in the 

30 art of progressive lens design from the geometrical properties of the surface without reference to the overall 
physical configuration of use of the lens. The second way of computing the purely optical quantitative 
measures of lens performance may be called true optical measure which takes such ray obliquity effects 
into account. To compute true optical measures for an ophthalmic lens, the actual physical relationship of 
the lens to the eye in use must be simulated, so that correct ray angles are obtained. Ray tracing 

35 techniques may then be used to compute the true optical measures of performance. 

In practice, a "family" of lens designs is constructed through a mathematical design model which 
contains a number of adjustable parameters that determine the quantitative performance measures of a 
particular design. Examples of such design models are given in the Applicants' EP-A-01 32955 referred to 
above. 

40 In accordance with the present invention, the quality of a particular design can be evaluated by a merit 
function rule to determine its merit figure. Thus the merit function may be regarded as a function of the 
adjustable parameters of the design model. In this context "optimization" refers to any process for 
determining values of the adjustable parameters which extremize, i.e., maximize or minimize, the value of 
the merit function. An optimum design is one which is determined by such optimal values of the adjustable 

45 parameters. 

The present invention provides a method for manufacturing an optimized ophthalmic progressive lens 
which has a distance viewing zone, a near viewing zone, and an intermediate zone in which a gradual 
change in optical power from the distance zone to the near zone occurs. The lens to be optimized, which 
has at least one quantitative performance characteristic, has its surface represented by a mathematical base 
so surface function. The base surface function is combined with a selected optimization function that includes 
at least one adjustable parameter and that may have at least one boundary condition imposed on it along 
the boundary of a selected subregion of the lens surface. A merit figure of a new surface function that is the 
combination of the base and optimization functions, which is related to the at least one quantitative 
performance characteristic by a selected merit function, is optimized through the at least one adjustable 
parameter of the optimization function. Through this optimization, there is obtained an improved surface 
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function which represents the surface of the improved progressive lens. A lens surface responsive to the 
improved surface function is then formed. 

The subregion of the lens surface can be selected to lie within the intermediate zone, or to be 
coextensive with that zone. Also, the at least one boundary condition imposed on the optimization function 

5 along the subregion's boundary may be homogeneous and these conditions may be imposed on the 
optimization function and certain of its partial derivatives. 

The selected merit function can be quadratic in the at least one adjustable parameter of the 
optimization function and can be either positive or negative definite, with the merit figure being optimized 
either when it is maximized or minimized. Also, the quantitative performance characteristic of the progres- 

w sive lens which is related to the merit figure by the merit function can be at least one of astigmatism, 
orthoscopy, mean curvature and deviation from a target surface. In a preferred embodiment, the perfor- 
mance characteristic is either surface astigmatism or a combination of surface and oblique astigmatisms. 

Further, the method of the present invention can be applied iteratively by varying the at least one 
adjustable parameter of the optimization function, evaluating the merit figure of the resulting new surface 

75 function that is the combination of the base and optimization functions, and repeating those steps until an 
optimum value of the adjustable parameter is obtained. Alternatively, the method can be applied so as to 
yield the optimum value directly. 

The present invention also embraces a multifocal, progressive lens having a near viewing zone, a 
distance viewing zone and an intermediate zone in which there is a gradual change in optical power along 

20 an eye path from the distance zone to the near zone. The progressive lens has the surface of its 
intermediate zone defined by a final surface function which is a sum of an optimization function and a base 
surface function. The optimization function takes on non-zero values in the intermediate zone which 
correspond to deviations in the height of a lens surface from the height of the base surface function, and 
includes at least one adjustable parameter which has its value derived by extremizing a merit figure. The 

25 merit figure is related to mean curvature or to surface astigmatism by a merit function and the optimization 
function changes the mean curvature or the surface astigmatism, respectively, of the final surface function 
by at least 0.1 reciprocal meters as compared with the base surface function and as measured at at least 
one point in the intermediate zone. 

The present invention further embraces a method for manufacturing an optimized ophthalmic progres- 

30 sive lens in which a base surface function is selected for a surface of a progressive lens having distance, 
near and intermediate zones. A gradual change in optical power from the distance zone to the near zone 
occurs in the intermediate zone. A merit function, which relates a weighted sum of squares of surface 
astigmatism at a plurality of preselected points to a merit figure of the lens, is constructed and an 
optimization function, which has at least one adjustable parameter, is chosen. The base surface function is 

35 combined with the optimization function to form a new surface function having the at least one adjustable 
parameter. The merit function of the new surface function is extremized to obtain an improved surface 
function that represents an improved surface of the progressive lens, and a lens surface responsive to the 
improved surface function is then formed. 

In a further embodiment a bounded subregion of the lens surface is defined and on the boundary zero 

40 is imposed on the optimization function and its first three partial derivatives across the boundary. The 
subregion can lie in the intermediate zone or can be coextensive with the intermediate zone. The 
optimization function can have the form of a polynomial in a curvilinear coordinate system with the 
coefficients of the polynomial including the at least one adjustable parameter. 

Embodiments of the present invention will now be described in greater detail by way of example only, 

45 with reference to the accompanying drawings in which: 

Figure 1 is a view of a progressive lens surface having multiple regions and showing a Cartesian 
coordinate system used to define the lens surface; 

Figure 2 is a graph of the surface height of a base surface function of a starting lens and of the 
surface height function of an improved lens provided by the method of the present invention; 

so Figure 3a is a view of a lens surface showing one location of a subregion; 

Figure 3b is a view of a lens surface showing another location of a subregion having computation 
points shown within it; 

Figure 4 is a view of a lens surface showing curvilinear coordinate curves used to define the lens 
surface; 

55 Figure 5 is a conceptual flow chart of a design procedure in accordance with the present invention; 

Figure 6 is a contour plot of the surface measure of astigmatism of a starting lens surface to be 
optimized by the method of the present invention; and 
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Figure 7 is a contour plot of the surface measure of astigmatism of an improved lens surface 
optimized by the method of the present invention. 



General Framework 

Referring to Figure 1, a progressive lens surface 1 may be considered as divided into three zones: a 
distance zone 10, a near zone 12, and an intermediate zone 14. The lens surface 1 can be described 
mathematically by giving its height above a reference plane. If points in the reference plane are labelled 
10 according to a Cartesian coordinate system having axes X and Y as indicated by the arrows shown, the 
height of the lens surface can be written as a mathematical function of points (x,y) in the reference plane. 

In proceeding to develop the optimization method of the present invention, it can be observed that a 
lens surface height function can also be expressed as a combination of two functions, thus: 
Z f (x,y) = Z b (x,y) + Z h (x,y) (1) 
75 where Z b (x t y) is a base surface function, Z h (x,y) is an optimization function and Z ( (x,y) is an improved 
surface function. Z b (x,y) represents a starting lens surface which is to be improved by the method of the 
present invention. In a preferred embodiment the starting lens surface may be one obtained by employing 
the methods set forth in Applicants' above-referenced U.S. patent application. The optimization function Z h - 
(x,y) is created by the present method so that a new lens surface represented by the combination of the 
20 base surface function and the optimization function provides a lens having quantitative measures of optical 
performance which are "improved" over those of the lens represented by Z b (x,y). The new lens surface 
having the optimal performance is represented by the improved, of final, surface function Z f (x,y). 

The base surface function may be continuous and its first several partial derivatives may be continuous 
as well. Similarly, the optimized lens represented by the improved surface function Z f (x f y) may also satisfy 
25 continuity conditions. As a result, the optimization function Z h (x,y) also satisfies similar conditions. 

An illustration of the effect of the optimization function Z h (x,y) on the height of the lens surface 
represented by the base function Z b (x,y) is shown in Figure 2. For clarity, the values of the two surface 
functions are illustrated for the Y coordinate only, i.e. for a constant X coordinate. The improved surface 
function Z f (x,y) is shown by the solid line, the base function Z b (x,y) is shown by the dotted line and the 
30 optimization function Z h (x,y) is shown by the dashed line. It can be noted that the optimization function can 
take on both positive and negative values. 

in practice it may be convenient to define Z h (x,y) such that the function takes on non-zero values only 
within a predetermined subregion of the lens surface and is zero everywhere outside the subregion. The 
location and shape of the subregion can be determined arbitrarily, such as the cross-hatched region shown 
35 in Figure 3a. in a preferred embodiment the subregion corresponds generally with the intermediate zone of 
a progressive lens, as shown by the subregion R in Figure 3b. Alternatively, the optimization function may 
take on non-zero values in a subregion which encompasses areas both within and outside the lens. 

However, when the subregion is restricted to only a portion of the lens surface, the continuity conditions 
introduced above require that the optimization function Z h (x,y) satisfy homogeneous boundary conditions; 
40 Z h (x,y) and some of its derivatives must be identically zero along the boundary of the subregion. Such 
boundary conditions cause Z h (x,y) to smoothly approach the zero value that it was defined to have outside 
the subregion. At least the optimization function and its first partial derivative must be zero along the 
boundary so that there are no sharp edges in the surface of the lens. In an embodiment described in more 
detail below, the optimization function and its first three partial derivatives are set to zero on the boundary. 
45 In accordance with the present invention, for any particular form of the optimization function selected at 
least one adjustable parameter will be provided. In general, the optimization function can be any function of 
the coordinates which is determined by one or more constant coefficients; such forms include polynomials 
as well as Fourier series. Each form of the optimization function corresponds to a "family" of lenses as 
described above. The adjustable parameters, which are the scalar coefficients of the optimization function, 
so also affect the improved surface function Z f (x,y). In the description below, the adjustable parameters are 
defined as the elements of an array A h . Accordingly, for this description, a more explicit functional notation 
is adopted for Equation 1: 

Z,(x,y; A h ) = Z b (x,y) + Z h (x,y; A h ) (2) 

It will be understood that the new surface function formed by the linear combination of Z b {x,y) and Z h - 
55 (x,y) shown in Equation 2 can also be formed by other, e.g. multiplicative, combinations of Z b and Z h . 

One way of representing the optimization function is by an expansion in terms of a system of functions 
of the form Pij(x,y), or more generally, Pj^u.v) where u and v are curvilinear coordinates. In a preferred 
embodiment the coordinates (u,v) may correspond to those coordinates (u,v) defined in Applicants' above- 
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referenced U.S. patent application. A point (u,v) in the curvilinear system can be expressed as a point (x,y) 
in the Cartesian system through a transformation: 

u = F(x.y) (3a) 

v = G(x.y) (3b) 

5 It will be understood that by using the coordinate transformation of Equations 3a and 3b, functional 
expressions in the Cartesian system are expressible in a different, and often more convenient form, in the 
curvilinear system. Figure 4 shows a view of a surface with Cartesian coordinate axes X and Y indicated by 
the arrows shown and an example of curvilinear coordinate curves 26 for constant v and coordinate curves 
28 for constant u. 

10 In terms of curvilinear coordinates the optimization function may be written as an expansion in terms of 
a system of functions Pij(u,v), thus: 

Z h (u,v;A h ) = J2 A h Pij(u,v) (4) 
ij 13 
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where the several A hi j are expansion coefficients, some of which can be zero. If the system of functions Py- 
(u.v) is a system of power functions of u and v, then Z h (u,v;A h ) is just a polynomial in u and v. The several 
A hiJ comprise the elements of the array A h of potentially adjustable parameters referred to in Equation (2). 
For the optimization function Z h (x,y; A h ) or its transformation Z h (u t v; A h ) to satisfy the special homogeneous 
boundary conditions, it is sufficient that each member of the system of functions Pjj(x,y) or Pjj(u,v) satisfy 
those conditions. 

Considering Equations 2 and 4, it can be observed that the improved lens represented by the improved 
surface function Z f (x,y; A h ), and therefore the quantitative performance measures of the improved lens, are 
effectively determined by the array of adjustable parameters A h . In accordance with the present invention, 
the adjustable parameters are varied so as to extremize a merit figure produced by a selected merit 
function. 

A merit function M(A h ) for a progressive lens may be generally defined as follows: 
M(A h ) = / R m(x,y.a,jS. A h ) dR (5) 
where M(A h ) provides a global measure of lens performance in a surface region R of the lens surface based 
on one or more quantitative performance measures, a,0,... evaluated at points (x,y) in the region R. 
Evaluation of the merit function M(A h ) for particular values of the array of adjustable parameters A h yields 
the merit figure for the lens defined by those parameters. The point merit function m(x,y,a,0, Ah) can be 
thought of as a local measure of lens performance at a point (x,y) on the lens surface. The point merit 
function m(x,y I ...;A h ) is the combination of quantitative measures of lens performance referred to above. 
Alternatively, the global merit function may be defined as a summation over a finite set of points in the 
region R (e.g. the small crosses in Fig. 3b), thus: 

M ( A h) = X>( x k/Yk' « / £ i ? A h ) (6) 

k 



Equation 6 is more amenable to numerical computation by computer. 

4 s in terms of merit functions such as those given in equations (5) and (6), "optimization" refers to any 
process for finding those values of the array A h of adjustable parameters which extremize the value of the 
merit function. Depending on the form of the merit function selected, optimization can involve either 
minimizing or maximizing the value of the merit function. As mentioned above, a variety of such methods 
are known to those skilled in the art. 

50 In accordance with one embodiment of the present invention, the values of the quantitative optical 
performance measures of astigmatism, mean curvature and orthoscopy at any given point are computed by 
surface measure, i.e., from the local properties of the lens surface as defined by the adjustable parameters, 
for use in the merit function. The relevant local properties of the lens surface are second derivatives of the 
surface function computed relative to a plane tangent to the surface at the point. These tangential partial 

55 derivatives can be denoted as Z^ t Zyy, and Z^ where Z** is the second tangential partial derivative with 
respect to x of a surface function Z(x,y), etc. Using this notation, the quantitative optical performance 
measures may be computed as follows: 
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Astigmatism: 
Mean Curvature: 

5 

Orthoscopy: 

It will be understood that by astigmatism is meant surface astigmatism which is the difference in the 
principal curvatures of the surface at the point of interest. Surface astigmatism as defined by Equation 7a 

70 can have units of reciprocal meters which can be related to the common ophthalmic units of diopters by a 
multiplicative factor of the index of refraction reduced by unity. Similarly, the mean curvature given by 
Equation 7b can also have units of reciprocal meters which are convertible to diopters of spherical 
equivalent power by the same multiplicative factor. Also, orthoscopy relates to the rectangularity of images 
which are observed through the lens. 

75 In one embodiment, the quantitative measures of lens performance are surface measures, described by 
Equations 7, which depend only on the geometrical properties of the lens surface. In another embodiment, 
the present method can employ true measures which are somewhat more difficult to compute. As 
discussed above, true optical measures such as oblique astigmatism, which arises when light rays 
encounter a refracting surface at angles off the normal to the surface, can be computed for the points in the 

20 optimization subregion by well-known ray tracing techniques. Normally, oblique astigmatism within the 
central thirty to forty millimeters of an ophthalmic lens is negligible, allowing the surface astigmatism alone 
to be used as a quantitative measure of lens performance. However, for points farther from the center of the 
lens, oblique astigmatism becomes more significant and so a net astigmatism, i.e. a combination of surface 
and oblique astigmatisms, can be used as a quantitative measure of lens performance in the present 

25 method. Of course it will be understood that net astigmatism can be used for all points in the optimization 
subregion, not merely for points away from the lens center. 

In accordance with the present invention, a class of measures called deviation-of-performance mea- 
sures permit the introduction of the desirable features of a known or target lens design into a new design. 
Differences between the surface measures (e.g.,a, 0, 7) of a lens design and corresponding surface 

30 measures of a target design are incorporated into the merit function in such a way that the value of the 
merit function, i.e. the merit figure of the improved design, decreases as the differences between the 
measured performances of the two designs decrease. In this way, for a merit function which is optimized by 
minimizing its value, the method of the present invention produces an improved lens design which 
approaches the optical performance of the target design. 

35 In accordance with the present invention, each of the quantitative deviation-of-performance measures is 
related to but not restricted to one of the optical measures of astigmatism, mean curvature and orthoscopy 
and the surface height of the lens. These relationships may be given by: 
5 = a - a t (8a) 
€ - 0 - jS, (8b) 

40 f = y - 7t (8c) 
U =* 2 f - Z t (8d) 

where 5 is the deviation-of-performance measure for astigmatism, a is the surface measure of astigmatism 
for the lens being optimized, and <* t is the surface measure of astigmatism of the target design. The 
deviation-of-performance measures € for mean curvature, f for orthoscopy and 1? for surface height are 
45 similarly defined. 

To emphasize that the quantitative performance measures depend on the particular design being 
considered (as specified by the array A h ) and on the particular point at which they are evaluated (as 
specified by its coordinates), the following functional expressions may be written: 



a 




<* (x.y; A„) 


(9a) 






£ (x.y; A h ) 


(9b) 
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y (x.y; A h ) 


(9c) 
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5 (x,y; Ah) 


(9d) 


€ 




€ (x,y; A h ) 


(9e) 






r (x,y; A h ) 


(9f) 






n (x.y; A h ) 


<9g) 



The merit function M(A h ) given by Equation 6 is a function of the quantitative performance measures 
given by Equations 9a-9g; thus combining these equations gives: 
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K(A h ) - *(*y:>YK> d , $ f Y . S ' <£ < Z >$> A n> do) 

which, as in Equation 5, could be written as an integral over the region R. 

A conceptual flowchart of the method of the present invention is illustrated in Figure 5. It can be noted 
that the particular sequence of the steps shown in the figure is substantially arbitrary and should not be 
considered the only acceptable sequence. A base surface function is selected which represents the surface 
of a progressive lens, the performance of which is to be optimized. It will be understood that the base 
surface function can represent a lens in accordance with the Applicants' above-cited U.S. patent application, 
or even a lens previously optimized by the present method. A subregion of the surface is selected for 
optimization and a convenient coordinate system is chosen. Typically, a generalized curvilinear coordinate 
system is chosen so that one or the other of the coordinates will be constant curves important in the 
optimization method, such as a curve defining the boundary of the optimization subregion or a curve 
defining an eye path from the distance zone to the near zone. The functional form pi the optimization 
function is then selected; the form can be any two-coordinate function which is determined by one or more 
scalar coefficients. Useful forms of the optimization function are polynomials, which are linear combinations 
of power functions of the coordinates, as well as Fourier and other types of series. Since the final surface of 
the optimized lens will be represented by a combination of the base surface function and the optimization 
function, the final surface is related to the coefficients of the optimization function selected. The coefficients 
are treated as adjustable parameters in the optimization method and their values are selected to optimize 
the performance of the lens having a surface represented by the combination of the base and optimization 
functions. 

The present method optimizes lens performance which is expressed in quantifiable performance 
characteristics, such as astigmatism, mean curvature, orthoscopy and deviation from a target surface. One 
or more of these performance characteristics, which generally have different values at each point on the 
lens surface, is selected and related by a merit function to a global merit figure for the lens. The particular 
form of the merit function is usually selected to be positive or negative definite to assure the existence of an 
extremum. Nevertheless, it is possible that indefinite merit functions can be used successfully by a skillful 
designer. 

After the form of the merit function is selected, the points in the optimization region where the 
performance of the lens is to be evaluated are chosen. The number and locations of these points are 
influenced by the computation time permitted and by the skill of the designer. Besides the number and 
location of points, weighting factors which control the importance of each of the performance characteristics 
are selected for each point. In this way the designer has direct control over the performance of the lens and 
can preserve in selected areas the performance of the starting lens if so desired. 

Depending on the form of the merit function selected, the values of the scalar coefficients of the 
optimization function which extremize the merit figure of the lens can then be determined either directly or 
iteratively. If the merit function has a form which has an extremum that is readily calculable from its 
coefficients, the optimization can proceed directly; an example of such a merit function is a quadratic form 
in the adjustable parameters. On the other hand, more complicated forms of merit functions can be used by 
applying the optimization iteratively, i.e. by selecting starting values of the adjustable parameters, calculat- 
ing the merit figure, varying the adjustable parameters, recalculating the merit figure, etc. until an extreme 
value (and the corresponding optimal values of the adjustable parameters) are obtained. 

Once optimal values of the adjustable parameters have been obtained for a given selection of 
optimization region, performance characteristics, points, weights, and merit function, the designer evaluates 
the resulting design and is free to vary these selections as indicated by the return arrows in Figure 5. 
Typically, the designer will more often vary the selection of points and weights, as indicated by the solid- 
line return arrow, than the other selections, as indicated by the dashed-line return arrows. The design 
process illustrated by the figure can continue until the designer obtains a satisfactory lens design. 

To illustrate how an optimization can proceed directly, consider a scalar quadratic function MF of a 
vector variable A: 

MF(A) = A T MA + B T A + C (I) 
where M is a given non-singular square matrix, B is a given vector and C is a given scalar. Without loss of 
generality, M may be taken as symmetric, but the matrix M must be definite, i.e. the term A T MA must be 
either non-negative for all values of A (positive definite), or non-positive for all values of A (negative 
definite). This requirement assures the existence of an extremum. 
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The value of A which extremizes the function MF is to be determined. For now, it will be assumed that 
M, B and C are known. In order to find an extremum of MF, its gradient with respect to the independent 
variable vector A is taken: 

GRAD MF = 2MA + B (II) 
Setting the gradient equal to zero and solving for the extremizing value of A gives: 

A mln = -(1/2)M~ 1 B (111) 
In practice, the values of M, B and C are unknown, but how to compute MF(A) for any value of A is known. 
M, B and C will be determined by evaluating MF on a suitably selected set of vectors H. Once M, B, and C 
have been found, equation (III) may be used to find the extremum of MF. 

A set of vectors H where Hi is the vector with exactly one non-zero component in the i'th position can 
be defined by: 

H, = (0,0, .... hu 0, 0) , i = 1 ...n (IV) 

where hj is the value of the non-zero component and n is the number of components of A. 

To compute M, B and C, it is sufficient to compute their components. First evaluate MF on the zero 
vector; this gives the value of C: 

MF(0) = 0 T MO = B T 0 + C = C (V) 
Next, evaluate MF on Hi: 

MF(H|) = Hj MHj + B T H| + C = h? M„ + h,Bi + C (VI) 
and then on (Hj + Hj): 

MF(Hi + Hj) = (Hi + Hj) T M(Hj + Hj) + B T (H, + Hj) + C = hf M n + 2hihjMij + h?M ii + h,Bi + hjBj + 
C (VII) 

Using Equations (V) and (VI), this can be written: 

MF(Hj + Hj) = MF(Hj) + MF(Hj) - MF(0) + 2h i H j M ii (VIII) 
Solving Equation (VIII) for M^ gives: 

Mij = [MF(Hj + Hj) - MF(H|) - MF(Hj) + MF(0)]/2hjhj (IX) 
This expression also gives M n . Thus, solving Equation (VI) for B t : 

Bj = [MF(Hi) - MF(0) - hf M n ] / h, (X) 
This completes the computation of the components of M, B and C from evaluations of MF on 0 and the set 
H. Equation (III) may now be used to compute the extremizing value of A. 



Specific Embodiments 

Having set out the general framework of the optimization method of the present invention, specific 
examples of a particular lens will be described. Ophthalmic progressive lenses have distance, intermediate 
and near viewing zones, in this example, the optimization method will be applied to optimize the 
performance of a progressive lens in the intermediate zone substantially as shown in Figure 3b. 

If the optimization subregion R is limited to less than all the lens surface, it has already been disclosed 
that the optimization function Z h (x,y) must satisfy the homogeneous boundary conditions. In carrying out the 
calculations of the present method, it may be convenient to employ a coordinate transformation from 
Cartesian coordinates to curvilinear coordinates given by: 

u = F(x,y) = x (11a) 

v = G(x,y) = [f d (x)-y]/[f d (x)-f n (x)] (1 1 b) 
where f d (x) is a function, the projection of which onto the lens surface defines the top boundary 16 of the 
optimization subregion R, and f n (x) is also a function, the projection of which defines the bottom boundary 
18 of the subregion R. The functions f d (x) and f n (x) are arbitrary, but can be low-degree polynomials. In a 
numerical example described in more detail below, second-degree polynomials are used. With the 
coordinate transformation of Equations 11a and 11b, the optimization subregion R conveniently comprises 
those points (u.v) on the lens surface having a v-coordinate between zero and unity. 

As already noted in regard to Equation 4, in a curvilinear coordinate system the optimization function 
can be expressed as an expansion in a system of functions Pij(u f v) provided the members of the system 
satisfy the homogeneous boundary conditions. The members of a system of functions which meet those 
conditions are given by: 
Pij(u,v) = S(u,v)v i d (1-v) t„ ukO-S-vy (12) 

where S(u,v) is a multiplicative auxiliary function such as that disclosed in Applicants' above-cited U.S. 
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patent application, the disclosure of which has been incorporated by reference. Also, l d and l n are fixed 
exponents which provide that the system of functions Pij(u.v) satisfies the homogeneous boundary 
conditions, and i and j are integers. Using the system of functions having members defined by Equation 12, 
the optimization function Z h (u,v;A h ) given by Equation 4 can be written as: 

Z h (u,v;A h )=S(u,v)v ^ d (l-v) ^ n Z A h uico.S-v)^ (13) 

ij ij 

The multiplicative auxiliary function S(u,v) can take any form so long as the members of the system of 
functions Pij(u.v) satisfy the homogeneous boundary conditions. Recalling that the coordinate transformation 
of Equations 11a and 11b produces an optimization subregion having boundaries of v = o and v = 1, the 
members of the system of functions Pij(u,v) will satisfy the homogeneous boundary conditions so long as S- 
(u.v) does not include factors of v~ td or (1-v)~ tn . This is equivalent to a requirement that S(u,v) be finite on 
the boundary. The exponents d and in provide that the partial derivatives of Pij(u,v) of orders less than id 
and in exist and are zero on the boundary. Suitable selection of the multiplicative auxiliary function 
provides a degree of control to the skillful lens designer. 

In the method of an embodiment of the present invention, optimization of the performance of a lens is 
accomplished by varying one or more adjustable parameters that determine the shape of the lens in such a 
way that the value of a merit function is minimized. Alternatively, the merit function can be so constructed 
as to indicate optimized lens performance when its value is maximized. In accord with another embodiment 
of the present invention, the point merit function m(x,y,a,/3,'yA<?,iM?:Ah) is preferably constructed to be 
definite to help ensure the existence of a minimum or maximum in the merit figure. By definite it is meant 
that the point merit function is positive definite or negative definite. If positive definite, it takes on only non- 
negative values; if negative definite it takes on only non-positive values. 

One well-known kind of positive definite function which can be used as the point merit function is a sum 
of squares in which the relative importance of each element of the sum is determined by a positive 
multiplicative weighting factor. In that case, the point merit function can be written as: 

m(x,y l o,i8,7,5 t € f ?,i ? ;A h ) = Wa a2 + w ) 3 + w y + ws 62 + w € f2 + w^ { ' 2 + w^ *> 2 (14) 
where the w f s are the weighting factors. In general, the values of weighting factors can be dependent on the 
point (x,y) on the lens surface where the point function is to be evaluated. It will be appreciated that in this 
way the relative importance of the quantitative measures of lens performance at each point in the 
optimization subregion R can be selectively varied and adjusted by the designer. Thus, the designer is 
provided with direct control over selected performance characteristics, especially in selected regions of the 
lens. 

In another embodiment, it is possible to define as zero in Equation 14 the weighting factors of ail the 
quantitative performance measures other than astigmatism «. Such a definition results in a point merit 
function and a global merit function which depend only on the surface measure of astigmatism; thus, from 
Equations 9, 10 and 14: 

M(A h ) = / R m(x,y,a;A h )dR = / R w a . a2 (x,y;A h )dR (15) 
or equivalently: 
M(A h ) = 

X 

o2 (x k ,y k ;A h ) (16) 

where the positional dependence of the weighting factor w a has been made explicit and the index k 
identifies a predetermined set of points (xi.yt). (x2,y 2 ),...,(xK,yK) in the optimization subregion R. These 
points 22 are shown in Figure 3b by the crosses in the intermediate region 14. 
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Numerical Example 

The following description illustrates the operation of the method of the present invention by a specific 
numerical example. Numerical values of the improved surface function Zf(x,y;A h ), given by Equation 2, 

5 which were generated by the present method for each point (x,y) are tabulated in Table 1 of the Appendix. 
Also given in Table 1 are values of the surface measure of astigmatism a (x,y) which are computed from 
Equation 7a for each point of the improved surface. An azimuthal angle, Q (x,y), of the surface astigmatism 
which gives the angle with respect to a reference direction of the larger of the principal curvatures is also 
listed. Similarly, numerical values of the starting surface represented by the base function Z b (x,y) are 

10 tabulated for each point in Table 2 of the Appendix. 

In the example, the optimization function Z h (u,v;A h ) was given by Equation 13 with t d = l n = 4 in the 
optimization region R; thus, the optimization function and its first three partial derivatives were zero on the 
boundary. Outside the optimization region, the optimization function was identically zero. The form of the 
multiplicative auxiliary function S(u,v) used in the example was: 

15 S(x,y) = exp [-p(x,y)-q(x,y)] (1 7) 

where the auxiliary function is expressed in the Cartesian coordinate system for convenience. The 
component functions that comprise the auxiliary function used are given by: 

P(x,y) - [<y-h)/(h c -h)] *h (x/a*)* (1 8) 
and: q(x,y) = [<y-y c -y m r x)/ayp for x £ 0 (1 9a) 

20 = [(y-y c -y m l x)/oy ] 3 for x<0 (1 9b) 

The constants h, h c ,n hf a Xl y c ,y mrt cy and y m i in Equation 18, 19a and 19b are tabulated in Table 3 of the 
Appendix, and control the shape and strength of the multiplicative auxiliary function. In this way, the 
multiplicative auxiliary function provides to the skillful designer a degree of direct control in addition to that 
provided by the weighting factors in the point merit function. The designer can achieve a desired overall 

25 lens performance by shaping the result computed with the selected weighting factors by the selected 
multiplicative auxiliary function. 

The functions f d (x) and f n (x) used in the coordinate transformation of Equations 11a and 11b are given 

by: 

f d (x) = adX 2 + b d x + c d (20a) 
30 f n (x) = a n x 2 +b n x + c n (20b) 

which are second-degree polynomials. The coefficients of Equations 20a and 20b are listed in Table 4 of 
the Appendix. 

In this numerical example, a merit function of the form of Equation 16 is selected and evaluated in a 
predetermined set of K = 96 points. The Cartesian coordinates of the points and the weight assigned to each 
35 point are listed in the Appendix in Table 5. Also listed in Table 5 are the values of the surface measure of 
astigmatism a which are computed for the final surface Zf(x,y) according to Equation 7a. 

Finally, the numerical values of the array A h of adjustable parameters which are generated by the 
method of the present invention are listed in Table 6 of the Appendix. It can be noted that the optimization 
function produced was a sixth-degree polynomial in x and only a second-degree polynomial in y. The 
4Q degree of the optimization function is determined by a balance between the complexity of the surface 
needed to achieve a desired performance and the cost of carrying out the optimization. 

The improvement in the astigmatism performance of the lens which is achieved by the method of the 
present invention for the numerical example is shown by comparing Figure 6 and Figure 7. Figure 6 shows 
a contour plot of the surface measure of astigmatism a for the base function Z b (x,y) listed in Table 2. Figure 
45 7 shows a contour plot of the surface measure of astigmatism a for the improved surface Z f (x.y) listed in 
Table 1. Each contour represents a change in astigmatism of 0.5 diopter. The lesser magnitude of the 
astigmatism of the improved design over the starting design is evident from Figures 6 and 7. 

The preferred embodiments of the present invention can provide an improved method of progressive 
lens design and an improved progressive lens. 
so The preferred embodiments of the present ivention can also provide an efficient and interactive design 
procedure for ophthalmic progressive lenses, and, also, an optimization procedure for ophthalmic progres- 
sive lenses which is based on quantitative measures of lens performance. 

With the preferred embodiments of the present invention, the optimization procedure can permit the 
introduction of the desirable features of an existing lens design into a new design developed through the 
55 procedure. 

The preferred embodiments of the present invention can further provide a quantitative merit function 
relating to optical performance characteristics of a progressive lens. 

The preferred embodiments of the present invention can still further provide an optimization procedure 
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for progressive lenses which provides the designer with direct control over selected performance char- 
acteristics, especially in selected regions of the lens. 

It is, of course, possible to embody the invention in other specific forms than those of the preferred 
embodiment described above. This may be done without departing from the essence of the invention. The 
preferred embodiment is merely illustrative and should not be considered restrictive in any way. The scope 
of the invention is embodied in the appended ciaims rather than in the preceding description and all 
variations and changes which fall within the range of the ciaims are intended to be embraced therein. 
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Table 1: 


Improved 


Surface 


Function 




X 


v 


Zf fx, V) 


fx.v^ 


fx, 


-10.00 


30.00 


5.92 


0.08 


89.78 


- 5.00 


30.00 


5.45 


0.16 


94 . 42 


0.00 


30. 00 


5.30 


0.19 


88.35 


5.00 


30. 00 


5.45 


0.17 


80. 07 


10.00 


30. 00 


5.92 


0.10 


69.83 


-20.00 


25.00 


6.08 


0.02 


135. 05 


-15.00 


25. 00 


4.99 


0.17 


117. 06 


-10.00 


25. 00 


4.23 


0.28 


110.82 


- 5.00 


25/00 


. 3-77 


0.27 


100.82 


0.00 


25.00 


3.62 


0-26 


90. 00 


5.00 


25.00 


3.77 


0.27 


79 . 18 


10.00 


25. 00 


4.23 


0.28 


69.18 


15.00 


25.00 


4.99 


0.27 


60.28 


20.00 


25.00 


6.08 


0.16 


45.56 


-25.00 


20.00 


6.08 


0.06 


146.13 


-20.00 


20.00 


4.68 


0.20 


132.62 


-15.00 


20.00 


3 .62 


0.26 


125. 69 


-10.00 


20.00 


2.88 


* 0.23 


115.8 0 


- 5.00 


20. 00 


2.43 


0.21 


103 . 66 


0.00 


20.00 


2.29 


0.20 


90.00 


5.00. 


20.00 


2.43 


0.21 


76.34 


10.00 


20.00 


2.88 


0.23 


64.20 


15.00 


20. 00 


3 .62 


0.26 


54.31 


20.00 


20.00 


4.68 


0.29 


46.62 


25.00 


20.00 


6.08 


0.21 


39.53 


-25.00 


15.00 


4.99 


0.17 


146.33 


-20.00 


15.00 


3.62 


0.26 


141.94 


-15.00 


15.00 


2.58 


0.22 


134.14 


-10.00 


15.00 


1.85 


0.17 


123.12 


- 5.00 


15.00 


1.42 


0.14 


108.15 


0.00 


15.00 


1.27 


0.13 


90.00 


5.00 


15.00 


1.42 


0.14 


71.85 


10.00 


15.00 


1.85 


0.17 


56.88 


15.00 


15.00 


2.58 


0.22 


45.86 


20.00 


15.00 


3.62 


0.26 


38 . 06 


25.00 


15. 00 


4.99 


0.29 


32.52 


-30.00 


10. 00 


5.92 


0.17 


81.95 


-25.00 


10.00 


4.23 


0.06 


137.84 


-20.00 


10.00 


2.88 


0.18 


151.43 


-15.00 


10.00 


1.85 


0.17 


145.75 


-10.00 


10.00 


1.13 


0.12 


134. 63 


- 5.00 


10.00 


0.70 


0. 08 


116.38 


0.00 


10.00 


0.56 


0.06 


90.00 
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Table 1: 


Improved 


Surface 


Function 


z f (x f y) 


X 


Y 


Zf (x f y) 


(x r y) 


fx,v> 


-20.00 


-10.00 


3 .21 


1.10 


73.66 


-15-00 


-10.00 


2.13 


1.49 


73.11 


-10.00 


-10. 00 


1.34 


1.73 


72. 10 


- 5.00 


-10 o 00 


0.84 


1.43 


69.76 


0.00 


-10. 00 


0.65 


0.62 


58.58 


5.00 


-10. 00 


0.81 


0.65 


132.78 


10.00 


-10 . 00 


1.33 


1.81 


119-76 


15.00 


-10. 00 


2.16 


2.20 


115.50 


20,00 ' 


-10. 00 


3.28 


2.24 


115. 05 


25.00 


-10.00 


4.71 


2.17 


112.99 


30.00 


-10.00 


6.45 
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-25.00 
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5.61 


1.00 


87.59 


-20.00 


-15.00 


4.18 
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-15.00 


-15.00 


3,0*6 
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83.03 


-10.00 


-15.00 
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1.67 
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- 5.00 


-15. 00 


1.72 
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0.00 


-15.00 


1.51 


0.28 
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5.00 


-15.00 
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0.14 
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0.15 
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-10. 00 


-30.00 
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1.73 


100.74 


- 5.00 
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0.87 


109.57 
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Table 1: Improved Surface Function Z f (x,y) 



0.00 -30.00 6.58 0.26 143-47 

5.00 -30.00 6.76 0.33 1.09 

10.00 -30.00 7.35 0.16 31.77 
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Table 2: 


Starting 


Surface 


Function 


S3 v A / J: J 


X 


v 


Zft fx, v) 


fx, 


f X . 


-10.00 


30.00 


5.92 


0.15 


139 . 02 


- 5.00 


30.00 


5.45 


0.05 


105. 68 


0.00 


30.00 


5.30 


0.11 


80 . 22 


5.00 


30.00 


5.45 


0 .17 


71. 82 


10.00 


30. 00 


5.92 


0.22 


65 . 40 


-20.00 


25.00 


6.08 


0.29 


115 .39 


-15.00 


25,00 


4.99 


0.28 


120 . 48 


-10.00 


25.00 


4.23 


0.28 


110. 82 


- 5.00 


25.00 


3 .77 


0.27 


100. 82 


0.00 


25. 00 


3 . 62 


0.26 


90. 00 


5.00 


25.00 


3 .77 


0.27 


79 . 18 


10.00 


25.00 


4.23 


0.28 


69.18. 


15.00 


25.00 


4.99 


0.28 


60.54 


20.00 


25.00 


6.08 


0.14 


61.51 


-25.00 


20.00 


6.07 


0. 18 


96. 10 


-20.00 


20.00 


4.68 


0.24 


129 . 65 


-15.00 


20.00 


3.62 


0.26 


125 . 69 


-10.00 


20.00 


2.88 


0.23 


115 . 80 


- 5.00 


20,00 


2 .43 


0.21 


103 .66 


0.00 


20.00 


2.29 


* 0.20 


90.00 


5. 00 


20.00 


2 .43 


0.21 


76.34 


10.00 


20.00 


2.88 


0.23 


64.20 


15.00 


20.00 


3.62 


0.26 


54 .31 


20.00 


20c00 


4.68 


0.29 


46. 62 


25.00 


20.00 


6.08 


'0,14 


30.28 


-25.00 


15.00 


4.99 


0.12 


148.80 


-20.00 


15.00 


3.62 


0.26 


141. 94 


-15.00 


15.00 


2.58 


0.22 


134.14 


-10.00 


15.00 


1.85 


0-17 


123 . 12 


— 5.00 


15.00 


1.42 


0.14 


108.15 


0. 00 


15.00 


1.27 


0.13 


90. 00 


5.00 


15.00 


1.42 


0.14 


71. 85 


10.00 


15.00 


1.85 


0 . 17 


56 . 88 


15o00 


15.00 


2.58 


0-22 


45 .86 


20.00 


15.00 


3.62 


0.26 


38 o 06 


25.00 


15.00 


4o99 


0.29 


32 . 47 


-30.00 


10.00 


5.92 


0 o42 


14 .42 


-25.00 


10.00 


4.23 


0.20 


156 . 57 


-20.00 


10.00 


2.88 


0.21 


152 .08 


-15.00 


10.00 


1.85 


0.17 


145 . 75 


-10.00 


10.00 


1.13 


0 .12 


134 . 63 


- 5.00 


10.00 


0.70 


0.08 


116.38 


0.00 


10.00 


0.56 


0.06 


90. 00 


5.00 


10.00 


0.70 


0.08 


63.62 


10.00 


10.00 


1.13 


0.12 


45.37 


15.00 


10.00 


1.85 


0.17 


34.25 


20.00 


10,00 


2.88 


0.23 


27.34 


25.00 


10.00 


4.23 


0.16 


25.06 
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Table 2: 


Starting 


Surface 


Function 


Zb(*#Y> 


X 


Y 






fXrY) 


30 . 00 


10 . 00 


5.92 


0.20 


22.53 


-30 . 00 


5 . 00 


5.46 


0.37 


177.21 


-25 . 00 


5 . 00 


3 .78 


0.21 


111.88 


-20 . 00 


5 . 00 


2.44 


0.27 


90.26 


-15,00 


5.00 


1.42 


0.31 


81-95 


-10*00 


5 . 00 


0.70 


0.24 


81.21 


- 5.00 


5. 00 


0.28 


0.10 


85.43 


0 • 00 


5.00* 


0.14 


0.03 


85. 39 


5 . 00 


5 . 00 


0.28 


0.07 


91. 08 


10 . 00 


5 . 00 


0.70 


0.38 


101.24 


15.00 


5. 00 


1.42 


0.82 


102.33 


20 . 00 


5. 00 


2.44 


0.96 


98.98 


25 . 00 


5. 00 


3 .78 


0.86 


92.11 


30.00 


5.00 


5.47 


0.67 


81.67 


-30 . 00 


0. 00 


5.32 


0.51 


152.47 


-25 . 00 


0. 00 


3.65 


0.61 


109.55 


-20.00 


0 . 00 


2.32 


0.73 


95.79 


-15 . 00 


0.00 


1.30 


0.87 


84. 64 


-10.00 


0.00 


0.58 


0.95 


74.23 


- 5.00 


0.00 


0.15 


\ 0.84 


63.66 


0.00 


0. 00 


0.00 


0.20 


63.76 


5.00 


0.00 


0.15 


0.94 


128.81 


10 . 00 


0.00 


0.59 


1-55 


120. 68 


15.00 


0. 00 


1.32 


1.56 


109.33 


20 . 00 


0.00 


2.35 


1.50 


97.48 


25.00 


0.00 


3 .70 


1.51 


88.53 


30.00 


0 . 00 


5.38 


1.46 


81. 64 


-30.00 


- 5.00 


5.52 


0.91 


144.12 


-25.00 


- 5.00 


3.86 


0.85 


116.20 


-20.00 


- 5.00 


2.53 


0.92 


101.60 


-15.00 


- 5.00 


1.51 


1.09 


86.59 


-10.00 


- 5.00 


0.78 


1.39 


73.54 


- 5.00 


- 5.00 


0.32 


1.30 


65.74 


0.00 


- 5.00 


0.16 


0.36 


62.93 


5.00 


- 5.00 


0.32 


0.93 


135.45 


10.00 


- 5.00 


0.79 


1.75 


126.56 


15 - 00 


- 5.00 


1.57 


1.74 


119.41 


20.00 


- 5.00 


2.63 


1.45 


111.93 


25.00 


- 5.00 


3 .99 


1.34 


99.58 


30.00 


- 5.00 


5.67 


1.44 


88.68 


-30.00 


-10. 00 


6.04 


1.32 


140.44 


-25.00 


-10.00 


4.40 


1.07 


121.23 


-20.00 


-10. 00 


3 .08 


1.30 


107.41 


-15.00 


-10.00 


2.06 


1.58 


93-20 


-10.00 


-10.00 


1.31 


1.91 


84.55 


- 5.00 


-10.00 


0.84 


1. 68 


76.41 
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Table 2 : 


Starting 


Surface 


Function 


Zb( x 'Y) 


X 


Y 


. Z b fx r y) 


(x,y) 




0 . 00 


—10 . 00 


0-65 


0.66 


60 . 68 


5 . 00 


—10 .00 


0-81 


0-64 


133 .26 


10 . 00 


—10 • 00 


1-32 


1- 91 


117 . 12 


15 • 00 


—10 .00 


2-15 


2 - 20 


107 .58 


20- 00 


—10 - 00 


3 . 25 


1. 81 


105 .85 


25 o 00 


—10 - 00 


4-65 


1. 53 


107 .45 


3 0. 00 


—10 .00 


6.36 


1- 47 


105 . 68 


— 25 • 00 


—15 .00 


5 . 27 


1. 63 


120 .23 


—2 0 . 00 


—15.00 


3 . 97 


1. 98 


111.04 


—15 • 00 


—15 • 00 


2.96 


2 • 44 


102 . 05 


—10 • 00 


— 15 .00 


2-21 


2 - 69 - 


92.18 


— 5 « 00 


—15 - 00 


1. 72 


2 . 05 


74 . 95 


0 • 00 


—15 . 00 


1.51 


0.30 


59 .38 


5 o 00 


—15 . 00 


1 . 67 


0 - 13 


146 .57 


10 . 00 


—15 - 00 


2 . 22 


1- 49 


121.76 


15 .00 


—15 • 00 


3 . 09 


2.96 


106 .76 


20 . 00 


—15 . 00 


4 . 24 


2-31 


93 . 64 


25 . 00 


-15 o 00 


5-67 


1-76 


98 . 13 


—25 . 00 


— 20 . 00 


6 - 52 


2 • 47 


114 .50 


-20*00 


-20 - 00 


5-24 


2 . 77 


111. 06 


-15 . 00 


— 20 . 00 


4 . 25 


3.40 


106 . 21 


-10. 00 


—20 . 00 


3-51 


3-53 


86.18 


— 5.00 


-20 - 00 


2 . 99 


1- 33 


78 , 03 


0 - 00 


—20 - 00 . 


2 . 77 


0 . 13 


108 .25 


5 • 00 


—20 o 00 


2 . 93 


0 . 07 


13 - 60 


10 . 00 


—20 - 00 


3 . 49 


0. 36 


94 . 08 


15 • 00 


—20 - 00 


4 - 43 


3 . 23 


113 . 69 


20 . 00 


—20 - 00 


5 • 64 


3.20 


99 . 36 


25 . 00 


—20 - 00 


1 • 11 


2 . 09 


94 . 22 


—20 .00 


—25 . 00 


6 . 94 


3 - 47 


109 -20 


—15 o 00 


—25 - 00 


5 . 98 


4 . 03 


95 . 63 


—10 . 00 


—25 - 00 


5-23 


4 - 21 


80 . 74 


— 5-00 


-25 . 00 


4-67 


0 - 99 


102.97 


0 • 00 


—25 . 00 


4 - 45 


0 . 13 


116.57 


5 • 00 


-25-00 


4-62 


0-19 


9 .70 


10 o 00 


—25 -00 


5-19 


0 • 16 


82 . 89 


15 • 00 


—25 .00 


6 - 16 


1 • 18 


97 . 17 


20-00 


-25.00 


7.48 


4-97 


111.31 


-10.00 


-30.00 


7-38 


3-03 


95.03 


- 5-00 


-30.00 


6-81 


1.56 


96.13 


0-00 


-30.00 


6-59 


0-58 


90.79 


5-00 


-30.00 


6-76 


0-21 


69.94 


10.00 


-30.00 


7-35 


0.31 


79.00 
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Table 3 : Auxiliary Function Constants 
o- x = 45 -° C>~y = 25.0 

h = -50,0 y c = 5*0 

h c 8=8 0-0 Yml = °- 75 

n h = 2 Ymr = -0,80 



J5 



20 



30 



35 



Table 4: Coefficients of the Boundary 
Functions f d (x) and f n (x) 



25 a d = 0.011073 

b d = 0,000000 
c d = 5,000000 



a n = -0,064000 
b n = 0,360000 



C n = -19,000000 



40 



45 



50 
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Table 5: Point Merit Function m(x,y, JA^) 



Point 



Index 






Weight 


cc 


k 






W^fXp^Vy^ 


1 


10,00 


-25.00 


0.0000 




2 


10-00 


-19.00 


0.3600 


1.30 


3 


10.00 


-13 . 00 


0.5200 


2.02 


4 


10.00 


- 7.00 


0.6800 


2.00 


5 


10.00 


- 1.00 


0.8400 


1.61 


6 


10.00 


5.00 


0. 0000 


- 


7 


12 . 5 0 


-25.00 


0. 0000 




8 


12 .50 


-19.00 


0.4000 


2.33 


9 


12.50 


-13.00 


0.6500 


2.16 


10 


12.50 


- 7.00 


0.5926 


2.10 


11 


12.50 


- 1.00 


0.7320 


1.64 


12 


12.50 


5.00 


0.0000 




13 


15.00 


-25.00 


0.0000 




14 


15.00 


-19.00 


0.4000 


2.31 


15 


15.00 


-13.00 


0.0500 


2.20 


16 


15.00 


- 7.00 


0.5051 


2.04 


17 


15.00 


- 1.00 


0.6240 


1. 64 


18 


15.00 


5.00 


0.7429 


0.51 


19 


17.50 


-25. 00 


0. 0000 




20 


17.50 


-19.00 


0.2000 


1.97 


21 


17.50 


-13.00 


0.1500 


2.26 


22 


17.50 


- 7.00 


0^0000 




23 


17o50 


- lo 00 


0.5160 


1.66 


24 


17.50 


5.00 


0. 6143 


0.64 


25 


20.00 


-25.00 


0. 0000 




26 


20.00 


-19.00 


0. 0000 




27 


20.00 


-13.00 


0.0000 




28 


20.00 


- 7.00 


0.0000 




29 


20.00 


- 1.00 


0.4080 


1.69 


30 


20.00 


5.00 


0.4857 


0.69 


31 


22 .50 


-25.00 


0.0000 




32 


22.50 


-19.00 


0.0000 




33 


22.50 


-13.00 


0. 0000 




34 


22 .50 


- 7.00 


OcOOOO 




35 


22.50 


- 1.00 


0.3000 


1.75 


36 


22 .50 


5.00 


0.3571 


0.69 


37 


26.00 


-17.00 


0.1500 


2.26 


38 


25.00 


-19 .00 


0. 0000 




39 


25.00 


-13.00 


0.0000 




40 


25.00 


- 7.00 


0.0000 




41 


25.00 


- 1.00 


0.1920 


1.81 


42 


25.00 


5.00 


0.2286 


0.67 


43 


27.50 


-25.00 


0. 0000 




44 


27.50 


-19.00 


0. 0000 




45 


27.50 


-11.00 


b.oooo 




46 


30.00 


- 7.00 


0.2000 


2.33 
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ruiiction m^x^y, o< 


' A h> 


Point 










Index 






" tixy xi L- 




k 







x k-^k-i 




47 


27.50 


- 1.00 


0 Oft 4 0 


X . O / 


48 


27.50 


5.00 


o i onn 

U . X u u u 


n a -5 

u • O J 


49 


-10 . 00 


—25 00 


n i 5nn 

U • iJ u u 




50 


— 9.00 


— i q nn 


n a n n n 

U t4UUU 


1 • 62 


51 


—10 . 00 


—i t on 

JLJ « UU 




T £T "7 
X . O / 


52 


—10.00 


— 7 nn 


n as o o 

U . O © u u 


x • o4 


53 


—10 . 00 


— i nn 


U . O *t U U 


1 • 13 


54 


—10 . 00 


5 on 


n noon 
u « uu u u 




55 


—12 . 50 


— ? 5 nn 

— > • U U 


n nnnn 
u • u u u u 




56 


—12 o 50 


-iq on 


n ^ n nn 

U • JUUU 


i.ol 


57 


-12.50 


—i ^ nn 


n 7o nn 
u • / u u u 


1 . 72 


58 


-12 . 50 


— 7 nn 

/ • U U 


n ao n n 
u • o u u u 


X • OJ 


59 


-12 .50 


— 1 00 


n 7*3 2 n 


t on 
x • z u 


60 


-12 . 50 


5 00 


n nnnn 

U . UUUU 




61 


— 15 on 


—2 5 no 

6J< UU 


U • UUUU 




62 


— i ^ nn 


— t q n n 


U • UUUU 






— i 5 no 

X -? « u u 


1j * UU 


U • UUUU 




64 


•15 no 


— 7 on 
/ • u u 


U • jUj! 


1 . 50 


65 


— i 5 nn 


— i on 

X • U U 


n *5 /i n 


1.23 


66 


—15 .00 


5 no 

3 « U \J 


n nnnn 
u . u u U U 




67 


—17 .50 


— 25 no 

^ W • U U 


n nnnn 

u . UUUU 




68 


-17 . 50 


—i q no 


u . UUUU 




69 


—17 SO 


— i *i on 

1J « UU 


u « UUUU 




70 


•i. / . w yj 


— 7 no 

/ . UU 


n nnnn 
U . UUUU 




71 


-17 50 


— x • u u 


U . OXoU 


1 • 24 


72 


—17 .50 


5 no 

•J • u u 


n AT A - * 


U . 3/ 


73 


—20 OO 


—2 5 no 

« UU 


n nnnn 
U • UUUU 




74 


—20 00 


— i q nn 

1j • UU 


n nnnn 
u . UUUU 




75 


—20 . 00 


— i *^ nn 


n nnnn 

u • UUUU 




76 


—2 0 00 


— 7 no 
/ • u u 


n nnnn 
u • uu u u 




77 


—20 . 00 


— i nn 
x • u u 


n Anon 
u • 4 U o U 


1 . 25 


78 


—20 . 00 


5 nn 

^ « u u 


U ■ ft O 3 / 


0.40 


79 


—22 . 50 


—25 nn 


n nnnn 

u • UUUU 




80 


—22 50 


—i q nn 

X-7 « U U 


n nnnn 
u . u u u u 




81 


—22 .50 


-13 no 

1J » UVJ 


n nnnn 

u • UUUU 




82 


—22 50 


— 7 no 

/ • u u 


n nnnn 
u . UUUU 




83 


—22 .50 


— i nn 

X . u u 


n nnnn 

u • UUUU 




84 


—22 . 50 


5 nn 

-J . u u 


n *5 cr7 1 


n c t\ 


85 


-25.00 


-25 . 00 


0 0000 

w . U v U U 




^86 


-25.00 


-19.00 


0.1200 


0.74 


87 


-25.00 


-13.00 


0.0000 




88 


-25.00 


- 7.00 


0.0000 




89 


-25.00 


- 1.00 


0.0000 




90 


-25.00 


5.00 


0.2286 


0.52 
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Table 5: Point Merit Function m(x,y, ck ;Aj 1 ) 



Point 
Index 



10 



15 



20 



Weight 



91 


-27.50 


-25.00 


0.0000 




92 


-27.50 


-19.00 


0.0000 




93 


-30.00 


- 7.00 


0.2500 


1.41 


94 


-27.50 


- 7.00 


0.0000 




95 


-27.50 


- 1.00 


0.0000 




96 


-27.50 


5.00 


0.1000 


0.53 



Table 6 



25 



Optimized Values of the Adjustable Parameter 

Matrix A h 



30 



35 



0 

1 

2 



-1.75X10 



-3 



2 . 17xl0~ 4 -2 . 3 1X10~ 6 



3 .20x10 



-8 



3.63x10 3 -7.18X10 4 



3 .09x10 



-6 



1.41x10 



-7 



-1.53X10*" 2 1.87xl0~ 3 -1.50xl0~ 5 -6.45xl0~ 7 
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45 



50 
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Claims 

1. A method for manufacturing an optimized ophthalmic progressive lens having a distance viewing 
zone, a near viewing zone and an intermediate zone in which there occurs a gradual change in optical 
power from the distance zone to the near zone, the ophthalmic progressive lens having at least one 
quantitative performance characteristic, comprising the steps of: 

representing a surface of the progressive lens to be optimized by a base surface function; 

selecting a merit function relating the at least one performance characteristic to a merit figure of the 

progressive lens; 

choosing an optimization function having at least one adjustable parameter; 

defining a subregion of the surface of the progressive lens, the subregion having a boundary; 

imposing at least one boundary condition on the optimization function along the subregion boundary; 

combining the base surface function and the optimization function to form a new surface function; 

optimizing the merit figure of the new surface function having the at least one adjustable parameter of 
the optimization function to obtain an improved surface function representing an improved surface of the 
progressive lens; and 

forming a lens surface responsive to said improved surface function. 

2. The method of claim 1, wherein the subregion lies, within the intermediate zone. 

3. The method of claim 1 , wherein the subregion is coextensive with the intermediate zone. 

4. The method of any one of claims 1 to 3, wherein the at least one boundary condition is homogenous. 
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5. The method of claim 4, wherein the at least one boundary condition includes the conditions that the 
optimization function and the first three partial derivatives of the optimization function across the boundary 
are zero on the boundary. 

6. The method of any foregoing claim, wherein the merit function is quadratic in the at least one 
5 adjustable parameter. 

7. The method of claim 6, wherein an optimum value of the at least one adjustable parameter is 
obtained by direct solution. 

8. The method of any one of claims 1 to 6. wherein an optimum value of the at least one adjustable 
parameter is obtained by an iterative procedure. 

70 9. The method of any foregoing claim, wherein the merit function is positive definite and the base 
surface function and the optimization function are combined through addition. 

10. The method of claim 9, wherein optimizing the merit figure comprises minimizing the merit figure. 

11. The method of claim 9, wherein optimizing the merit figure comprises maximizing the merit figure. 

12. The method of any foregoing claim, wherein the at least one performance characteristic is at least 
75 one of astigmatism, orthoscopy, mean curvature, and deviation from a target surface. 

13. The method of any one of claims 1 to 11, wherein the at least one performance characteristic is 
surface astigmatism. 

14. The method of claim 13, wherein a square of a surface astigmatism, a 2 , is used in evaluating the 
merit figure and wherein the square of the surface astigmatism is obtained by evaluating at at least one 

20 point the expression: 

" 2 = (Zyy - Z^) 2 + 4Z 2 xy 

wherein Z is the surface height of the new surface function, and the derivatives Zxy, Zyy and Z w are 
computed relative to a tangent plane at the point of evaluation. 

15. The method of any one of claims 1 to 11, wherein the at least one performance characteristic is a 
25 combination of surface astigmatism and oblique astigmatism. 

16. A multifocal, progressive lens with a lens surface having a near viewing zone, a distance viewing 
zone and an intermediate zone in which there is a gradual change in optical power along an eye path from 
the distance zone to the near zone; 

wherein the surface of the intermediate zone is defined by a final surface function which is a sum of an 
30 optimization function and a base surface function; 

wherein the base surface function provides a gradual change in optical power along an eye path from 
the distance zone to the near zone; 

wherein the optimization function has non-zero values in the intermediate zone which correspond to 
deviations in a surface height of the lens from a surface height of the base surface function; 
35 wherein the optimization function has at least one variable whose value is derived by extremizing a 

merit figure calculated from a merit function which merit function relates the merit figure to mean curvature 
or to surface astigmatism; and 

wherein the optimization function changes the mean curvature or the surface astigmatism, respectively, 
of the final surface function, as compared with the base surface function, by at least 0.1 reciprocal meters at 
40 at least one point in the intermediate zone. 

17. A method for manufacturing an optimized opthalmic progressive lens, comprising the steps of: 
selecting a base surface function for a surface of a progressive lens having a distance viewing zone, a 

near viewing zone and an intermediate zone in which there occurs a gradual change in optical power from 
the distance zone to the near zone; 
45 constructing a merit function relating a weighted sum of squares of surface astigmatism at a plurality of 
preselected points to a merit figure of the progressive lens; 

choosing an optimization function having at least one adjustable parameter; 

combining the base surface function and the optimization function to form a new surface function; 

extremizing the merit figure of the new surface function having the at least one adjustable parameter of 
so the optimization function to obtain an improved surface function representing an improved surface of the 
progressive lens; and 

forming a lens surface responsive to said improved surface function. 

1 8. The method of claim 1 7, further comprising the steps of: 

defining a subregion of the surface of the progressive lens, the subregion having a boundary; and 
55 imposing on the optimization function boundary conditions that the optimization function and the first 

three partial derivatives of the optimization function across the boundary are zero on the boundary. 

19. The method of claim 18, wherein the subregion lies in the intermediate zone. 

20. The method of claim 1 8, wherein the subregion is coextensive with the intermediate zone. 

23 
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21. The method of any one of claims 17 to 20. wherein the optimization function (Z h ) chosen is of the 
form 



2 h( u ' v ) = ij MjU 1 ^ 



wherein (u,v) are curvilinear coordinates, and wherein coefficients A Vi include the at least one adjustable 
parameter. 
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